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Abstract. The e — e , e — i, i — i and charge-charge static structure factors are calculated for alkali 
and Be 2+ plasmas using the method described by Gregori et al. in [I]. The dynamic structure factors 
for alkali plasmas are calculated using the method of moments [2], [3]- In both methods the screened 
Hellmann-Gurskii-Krasko potential, obtained on the basis of Bogolyubov's method, has been used taking 
into account not only the quantum-mechanical effects but also the ion structure [4]. 
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• Sh 1 Introduction 

^h ^. The structure and thermodynamic properties of alkali and 
f-^ 1 alkaline earth plasmas are of basic interest and of impor- 
i— i, tance for high-temperature technical applications. Near 
the critical point and at higher temperatures the mate- 
rials are in the thermodynamic state of strongly coupled 
plasmas. Here we will go far beyond the critical point to 
the region of nearly fully ionised plasmas which is T > 30 
kK for alkali and T > 100 kK for alkaline earth plasmas. 
The investigation of thermodynamic properties of alkali 
plasmas under extreme conditions is not only important 
for basic research. There are many applications, e.g. in ma- 
terial sciences, geophysics and astrophysics. Furthermore, 
these studies throw some light on the complex picture of 
phase transitions in metal vapors which play an outstand- 
ing role in technological applications. 

Over the recent decades a considerable amount of effort 
has been concentrated on the experimental [5] , [5] and the- 
oretical [7]-[H] investigation of the behavior of alkali metals 
in the liquid and plasma state expanded by heating toward 
the liquid-vapor critical point. High-temperatures alkali 
plasmas are widely applied in many technical projects. For 
instance, Li is an alkali metal of considerable technological 
interest. Lithium is planned to be used in inertial confine- 
ment fusion, solar power plants, electrochemical energy 
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Table 1. The ionization energies h (eV) of alkali and Be atoms 
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storage, magnetohydrodynamic power generators and in 
a lot of other applications. Recent advances in the field of 
extreme ultraviolet EUV lithography have revealed that 
laser-produced Li plasmas are source candidates for next- 
generation microelectronics |10| . For this reason we believe 
that the study of basic properties of alkali plasmas are of 
interest. In the previous work we studied Li + plasma [J]. 
In this work, we consider Li + , Na + , K + , Rb + , Cs + and 
Be 2+ plasmas. For simplicity of the calculations we take 
into account here only single ionization for alkali plasmas 
(n e = Ui) and doubled for beryllium plasma (n e = 2rij), 
where n e , n, are the concentrations of electrons and ions 
respectively. Li, Na, K, Rb, Cs atoms have one outer elec- 
tron and Be 2+ has two outer electrons. In the table [1] the 
ionization energies of alkali atoms and beryllium atom are 
presented. Correspondingly we will study temperatures 
around 30000.fr for alkali and 100000.Fr for Be 2+ plasmas, 
where most of outer electrons are ionized, but the rest core 
electrons are still tightly bound. 

Recently, X-ray scattering has proven to be a powerful 
technique in measuring densities, temperatures and charge 
states in warm dense matter regimes In inertial con- 
finement fusion and laboratory astrophysics experiments 
the system demonstrates a variety of plasma regimes and 
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of high interest are the highly coupled plasmas Fa > 1 
{r u = z 2 e 2 /{4e k B Tr u ) with r u = (S/Ann,) 1 / 3 being 
the average ion-ion distance, e is the electric elementary 
charge and z - the ionic charge) and the electron subsys- 
tem exhibiting partial degeneracy. Such regimes can be of- 
ten found during plasma-to-solid phase transitions. Recent 
experiments with a solid density Be plasma have shown 
a high level of ion-ion interaction and their interpretation 
must account for significant strong coupling effects. The 
present study is devoted to the study of the static (SSF) 
and dynamic (DSF) structure factors for alkali and Be 2+ 
plasmas at temperatures T > 30 kK and T > 100 kK 
respectively. The structure factor (SF) is the fundamen- 
tal quantity that describes the X-ray scattering plasma 
cross-section. Since the SF is related to the density fluctu- 
ations in the plasma, it directly enters into the expression 
for the total cross-section. In the case of a weakly cou- 
pled plasma SF can be obtained within the Debye-Huckel 
theory or the random phase approximation (RPA), while 
at moderate coupling the RPA fails to predict the correct 
spatial correlations. However, recent works by Gregori et 
al. [T] , [T2] and by Arkhipov et al. [T3] have shown that the 
technique developed in the classical work of Bogolyubov 
provides sufficiently reliable expressions of SF in moder- 
ately coupled plasmas. 

For the determination of static and dynamic structure 
factors one needs to have a screened pseudopotential as 
an essential input value. The semiclassical methods allow 
to include the quantum-mechanical effects by appropri- 
ate pseudopotentials which resolve the divergency problem 
at small distances. This method was pioneered by Kelbg, 
Dunn, Broyles, Deutsch and others [14] - [16]; later it was 
significantly improved in [T7]-[TH]. These models are valid 
for highly temperature plasmas when the ions are bare or 
there is no significant influence of the ion shell structure. 
In order to correctly describe alkali plasmas at moderate 
temperatures one needs to take into account the ion struc- 
ture. For example, to describe the behaviour of alkali plas- 
mas, the short range forces between the charged particles 
are of great importance. For alkali plasmas at small dis- 
tances between the particles deviations from the Coulomb 
law are observed which are mainly due to the influence 
of the core electrons. The method of model pseudopoten- 
tials describing the ion structure was pioneered by Hell- 
mann. Hellmann demonstrated, using the Thomas-Fermi 
model, that the Pauli exclusion principle for the valence 
electrons can be replaced by a non-classical repulsive po- 
tential |20j . This method was later rediscovered and fur- 
ther developed for metals by Heine, Abarenkov and Ani- 
malu [2T], [35]. Heine and Abarenkov proposed a model, 
where one considers two types of interaction: outside of the 
shell, where the interaction potential is purely Coulomb 
one and inside, where it is constant. Parameters of these 
model potentials were determined using the spectroscopic 
data. Later on different pseudopotential models were pro- 
posed. For a more detailed review we refer to [22]. All 
these models are characterized by one disadvantage. Their 
Fourier transforms (formfactors) are not sufficiently con- 
vergent when the Fourier space coordinate goes to infin- 



ity. Gurskii and Krasko [23] proposed a model potential 
which eliminates this problem and provides smoothness of 
the pseuopotential inside the shell giving its finite value at 
small distances. The first attempt to construct the model 
pseudopotential for alkali plasmas taking into account the 
ion structure was made in the works [23], [25], where the 
Hellmann type pseudopotentials were used. In this work 
we use Hellmann-Gurskii-Krasko pseudopotential model 
for electron-ion interactions and its modified version for 
ion-ion interactions. There is also a high interest to con- 
struct a pseudopotential model of particle interaction in 
dense plasmas taking into account not only the quantum- 
mechanical effects including the ion shell structure at short 
distances but also the screening field effects. In the work 
[3] the screened Hellmann-Gurskii-Krasko potential was 
derived using Bogolyubov's method as described e.g. in 
[25], [27] and [25]. 

We consider Li+, Na+, K+, Rb+, Cs+ and Be 2+ two- 
component (TCP) plasmas with the charges ze_ = — e + 
and masses m, >> m e and the densities n e — zn, (z = 
1,2). We will calculate here the TCP static and dynamic 
structure factors, including quantum effects and the ion 
shell structure using the Hellmann-Gurskii-Krasko pseu- 
dopotential (HGK). To determine the static and dynamic 
structure factors we use the screened Hellmann-Gurskii- 
Krasko potential obtained in j4|. The method which is 
used for the calculation of the static structure factor is 
the TCP hypernetted-chain approximation developed for 
the case of absence of the local thermodynamic equilib- 
rium (LTE) by Seuferling et al [2J5] and further discussed 
and extended for SSF by Gregori et al. [TJ while for the dy- 
namic - the method of moments applied to two-component 
plasmas [2J , [3J . We would like to underline again that the 
inclusion of both components into the theory and a cor- 
rect account of the short-range electron-ion interactions, 
is essential for the understanding of the structure factors 
in plasmas. 



2 Pseudopotentials taking into account the 
ion structure. Hellmann-Gurskii-Krasko 
potential 

Clearly, the simple Coulomb law is not applicable to de- 
scribe the forces between the charges in alkali plasmas 
since, at small distances, there are strong deviations from 
Coulomb's law due to the influence of core electrons. 

In many problems of atomic and molecular physics one 
can divide the electrons of the system into valence and 
core electrons. A number of important physical properties 
are determined by the valence electrons. In a series of pio- 
neering papers Hellmann attempted to develop a model in 
which the treatment of atoms and molecules in alkali plas- 
mas is reduced to the treatment of valence electrons [2U] . 
Hellmann demonstrated that the Pauli exclusion principle 
for the valence electrons can be replaced by a nonclassi- 
cal potential (repulsive potential) which is now called the 
pseudopotential. Hellmann's idea was to replace the re- 
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quirement of orthogonality of valence orbital to the core 
orbitals by the pseudopotential and it was this idea which 
made the respective mathematical calculations much sim- 
pler. 

For the actual purpose of atomic and molecular calcu- 
lations Hellmann suggested a simple analytic formula. Let 
if be the sum of electrostatic, exchange, and polarization 
potentials, representing the interaction between a valence 
electron and the core of an atom. Let tp p be the repulsive 
potential. The Hellmann potential ipn = <f~ + <f P may be 
expressed as: 



Aexp(— or), 



(1) 



Here z is the ionic charge of the core; that is, if the nu- 
cleus contains Z positive charges and the core contains N 
electrons then z = Z — N. The constants A and a are 
determined from the requirement that the potential ipu 
should reproduce the energy spectrum of the valence elec- 
tron as accurately as possible. Later several modifications 
were introduced by Schwarz, Bardsley etc. into the deter- 
mination of the Hellmann potential parameters without 
changing the basic analytic form of the potential [30] . For 
instance, Schwarz improved the determination of the po- 
tential parameters of Be + , Mg + , Ca + , Sr + , Zn + and Li, 
Na, K, Rb, Cu periodic families consistently obtaining a 
better fit to the empirical energy levels |31| . 

However, all above pseudopotentials have one draw- 
back. They are usually described in the r space by a dis- 
continuous function or have a relatively hard core as in 
the case of Hellmann's potential. As a result, their Fourier 
transforms (formfactors), f{k), at k — > oo, do not guaran- 
tee the convergence of series and integrals of the perturba- 
tion theory. Alternatively, Gurskii and Krasko constructed 
a pseudopotential model excluding this shortcoming by in- 
troducing a continuous in the r space pseudopotential. To 
include the smoothness of the obtained electron density 
distribution in an atom, Gurskii and Krasko proposed the 
following electron- ion model pseudopotential [23 , 32 : 
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■ exp 
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where Rc e i = r Cei r B and a are determined experimen- 
tally using the ionization potential and the formfactor of 
the screened pseudopotential at the first nodes of the re- 
ciprocal lattice. The parameter rc e i is defined as a cer- 
tain radius characterizing the size of the region of internal 
electron shells. If such measurements are not available, 
the second condition is replaced by the constraint that 
the pressure P = at zero temperature in the equilib- 
rium lattice. The magnitudes are given in the SI system 
of units. In this work values of a, rc e i are taken from |33j . 
Since the first two terms in @ are identical to the po- 
tential of Hellmann [20] , we call this pseudopotential the 
Hcllmann-Gurskii-Krasko potential. The results of calcu- 
lation of both the bound energy and the phonon spectra 



using the Hellmann-Gurskii-Krasko (HGK) potential were 
found in a good agreement with the experimental data and 
can be used in a wide range of investigation of thermody- 
namic properties of alkali plasmas. Unfortunately there 
are no available HGK parameters for the Be 2+ ion. That 
is why we looked for alternative e — i potentials with the 
determined for Be 2+ parameters. It is the Fiolhais et al. 
pseudopotential |34) : 



")=- 



1 



ATreor c R 



{l-(l + /3R)exp(-aR)} (3) 



that we used to describe Be 2+ plasmas. Here R = r/r C) r c 
being a core decay length, a > 0, f3 = (a 3 — 2a)/4(a 2 — 1) 
and A — a 2 /2 — a/3. In [33] there are two possible choices 
of parameters: the "universal" and the "individual". We 
made a fit of the "universal" parameters of HGK to the Fi- 
olhais et al. pseudopotential, which are a — 3.72, r = 0.22. 
In [34] the universal parameters were chosen to obtain the 
best agreement between calculated and measured struc- 
ture factors of alkali metals. In Fig. Q] the comparison be- 
tween the electron-ion Fiolhais et al., HGK and Coulomb 
potentials for Be 2+ plasma are shown. One can easily see 
that HGK almost coincided with the Fiolhais et al. po- 
tential. In Fig. [5^ we display the pseudopotentials <p e i{r) 
for different alkali plasmas and the HGK pseudopotentials 
tpu(r) are represented in Fig.[2p; notice that tp e i(r) possess 
a minimum. The Hellmann type pseudopotentials for al- 
kali plasmas were proposed in [33] , [35] and the Hellmann- 
Gurskii-Krasko model for the ion-ion interaction shown in 
Fig. [5p for alkali plasmas is the following: 



HGK 



(r) 



1 — exp — 



Rci 



47T£ R Cl 



exp 



Rci 



(4) 



The values of rcu, a are not given in literature, particu- 
larly, Tea is taken hypothetically as the doubled value of 
that taken for e — i interaction rcu — 2rc e i taking in this 
way both ions cores (closed shells) into account. We will 
study this in more detail and compare with the hard-core 
potential described in 25 . In Table [2] the parameters of 
the Hellmann-Gurskii-Krasko potential for alkali elements 
and the elements of the second periodic family are pre- 
sented. We note that ip HGK ei (r) potential describes the 
interaction of a valence electron with the corresponding 
ion core of a charge z and radius Rc e i, while tp HGK u (r) 
describes the interaction between two ion cores of a charge 
z with the same radius ren- 
in principle, the choice of the parameters for the ion- 
ion interaction should be based on methods similar to 
those in [33]. Furthermore, our calculations led us to the 
conclusion that the potential is not sensitive to the a pa- 
rameter of the ion-ion interaction. That is why a was taken 
the same as in the electron-ion potential. In Figures [3^,, 
b the comparison between the HGK, hard core (HC) [2"5] 
and soft core (SC) (eq. (5) in [35]) models is presented. 
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Fig. 2. Comparisons among the HGK potentials of different alkali plasmas (in atomic units). 

Table 2. The parameters of the Hellmann-Gurskii-Krasko potential in atomic units. Here, in the case z = 2 the given parameters 
correspond to the interaction of a double charged ion with an electron. 
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5.954 


3.362 


2.671 
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2.214 
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2.588 
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2.575 


2.870 


TCei 


0.365 


0.487 


0.689 


0.779 


0.871 


0.22 


0.427 


0.571 


0.644 


0.698 


rcu 


0.73 


0.974 


1.948 


1.558 


1.742 


0.44 


0.854 


1.142 


1.288 


1.396 



lowers is basically working with wave functions of a few 




particles and not the multiparticle wave functions of the 
solid state. For this reason we cannot see strong argu- 
ments against the two-particle interaction we employ in 
the plasma state. Of course this is a working assump- 
tion which needs further justification. Anyhow we are con- 
vinced that the application of pseudopotentials of Hcllmann- 
type is much nearer to reality than the use of pure Coulomb 
potentials or hard-core potentials. Further, we would like 
to argue that the experimental investigations of alkali met- 
als near the critical point did not evidence the existence of 
deep differences between the two particle interactions in 
the liquid and the gaseous state [5] , [5] . What is clearly dif- 
ferent are the multi-particle interactions, however multi- 
particle effects are less relevant at the densities we consider 



R=r/r here. 



Fig. 1. Comparisons among the e — i HGK, Fiolhais et al. and 
Coulomb potentials for Be 2+ (in atomic units). 

The pseudopotentials which are used in our calcula- 
tions were originally developed for applications in the elec- 
tronic band structure and binding energies in alkali met- 
als. However the derivation used by Hellmann and his fol- 



2.1 Screening of the Hellmann-Gurskii-Krasko 
potential 

Effective potentials simulating quantum effects of diffrac- 
tion and symmetry [14] - |16) as well as significantly im- 
proved potentials [12]- [19] are widely used to determine 
the thermodynamic and transport properties of semiclas- 
sical fully ionized plasmas. In particular, Deutsch and co- 
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workers [TS], Qj5] obtained the following effective interac- 
tion potential of charged particles in a plasma medium: 



'•pe 



exp 



A,. 



-fcs7Tn2exp 



ALln2 

/ . (5) 

where A ee = H/^/m e kBT is the electron thermal de-Broglie 
wavelength. In the work [28] it was proposed to use the 
e — e effective potential (|5|), along with the correspond- 
ing e — i and i — i potentials at short distances and the 
screened potential, treating three particle correlations, at 
large ones. The transition from one potential curve to an- 
other was realized at the intersection point by the spline- 
approximation method. 

The pseudopotential model was developed only 
for highly ionized plasmas. Since most experimental data 
available refer to partially ionized plasmas at moderate 
temperatures when the ions partially retain their inner 
shell, it is of a high interest to construct the pseudopoten- 
tial model which takes into account not only the quantum- 
mechanical and screening field effects but also the ion shell 
structure. In order to include the screening effects, in [4] 
we applied the method developed in [27] and [25J. In [25J 
the authors suggested the classical approach based on the 
chain of Bogolyubov equations |26) for the equilibrium 
distribution functions where the potential ([5]) was taken 
as a micropotential. In [3] the Fourier transforms of the 
screened HGK were derived, eq. (11)-(14) therein, using 
the e — i, i — i Hellmann-Gurskii-Krasko pseudopotentials 
(J2J), (0]) and e — e Dcutsch potential ([5]) as the micropo- 
tentials. 

The screened HGK potential was obtained in U in the 
following way. In the Fourier space the system of intcgro- 
diffcrential equations which determines the screened HGK 
potential turns into the system of linear algebraic equa- 



tions: 



&ab(k) = <Pab(k)- -;—= [n e ip ae (k)@ eb (k) + n 2 (p al (k)@ ib (k)] 
k B T 

(6) 

where a,b — i, e. Solving the system ([6]) for two-component 
plasma one can derive the following expressions for the 
Fourier transform ^ a b(fcl3of the pseudopotential ^ a b{f)- 



ze 2 (2a- l)R ce ik 2 - 1 
e A k 2 (l + k 2 R Cei 2 ) 2 ' 



(7) 



®ee(k) - ^ [ k 2 (l + k- 



k 4 r Di 



(l + fc 2 A„ 2 )(l+fc 2 .R cii 2 ) 2 



fc 2 A„ 2 ) 

(2a-l)i?. eei 2 fc 2 -l N ~ 
(l+k 2 R c „ 2 ) 2 



A 1 



(2a+l)Rcii k 2 +l 
k 2 r Dz 2 (l+k 2 R aii 2 ) 2 I eX P 



(-£)}, w 



+ 



FT 



(b\ - f (2a+l)fl c „ 2 fc 2 + l 

l W ~ s„ A \ k 2 (l + k 2 R ail 2 ) 2 

( (2a-l)R cri 2 k 2 - 
\ (l + fc 2 -R cel 2 ) 2 



(2a+l)fl cii 2 fc 2 + l 
(l+k 2 \ cc 2 )(l+k 2 R cU 2 ) 2 



/1 fc 2 r De 2 (l+*; 2 J?, ctl 2 ) 2 LX P ^ 4b J j 



(9) 



Here J"D e , ro» are the Debye radius of electrons and ions 
respectively, with 1/roj 2 — z 2 e 2 rii/ (£ £;_bT), 
l/r D 2 = e 2 n e /{e k B T), b = (A ee 2 In 2)" 1 , 



1 A correction is made with respect to a misprint in [4]: 
in eqs. (|8I10|) there should be ~ exp ^ — j^j instead of ~ 

cx p(-J&) 
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k B Tln2ir 3 / 2 b- 3 / 2 e Q /e 2 and 



A = 1 



1 , (2a+l)R aii 2 k 2 + l 

k 2 r Da 2 (l + k 2 \ ea 2 ) + k 2 r Dz 2 (l+k 2 R aii 2 ) 2 



(2a+l)R cU 2 k 2 + l 
{l + k 2 \ ee 2 ){l+k 2 R cU 2 ) 2 



^ 

roe' 



(2a+l)R €i 2 k 2 + l 
k 2 r Di 2 (l+k 2 R cU 2 ) 2 



( (2a-l)R cei 2 k 2 -l \ 
\ (1+k 2 R aci 2 ) 2 J 

*p(-£)- (io) 



the interaction between the particles (charged spheres) 
to account for ion structure, quantum diffraction effects 
i.e., the Pauli exclusion principle and symmetry. Quan- 
tum diffraction is represented by the thermal de Broglie 
wavelength A rs = H/^2p' rs k B T'j, s with p' rs = the 



The pseudopotential <5> a b{r) can be restored from (J7J- 
HOP by the Fourier transformation 



$ab{r) 



27T 2 7 



@ab(k)k sm(kr)dk 



(11) 



The present approximation is restricted to the constraint 
Fa < 1 due to the employment of the linearization pro- 
cedure in the derivation of the resolved system of integro- 
diffcrcntial equations. 

In order to compare with the alternative i — i pseu- 
dopotential models we considered the Dalgic potential. 
S.S. Dalgic ct al. determined the screened ion-ion poten- 
tial [36] on the basis of the second order pseudopoten- 
tial perturbation theory using the Fiolhais et al. potential 



*«(*) = 



4irsr)k 2 



where ip F ei (k) is the pseudopotential local form factor. 
Here, we use instead of the Fiolhais potential the HGK 
®u GK (k) potential with the fitted to those of Fiolhais et 
al. potential parameters. Notice that x(&) is the response 
of the electron gas, where the Lindhard response of a non- 
interacting degenerated electron gas x (&) an d the local 
field correction (LFC) G(k) enter the expression in the 
standard way; the LFC accounts for the interactions be- 
tween the electrons. We used the LFC which satisfies the 
compressibility sum rule at finite temperatures obtained 
by Gregori et al. for the strong coupling regime [12]. Sim- 
ilar calculations have also been carried out by E. Apfcl- 
baum who used the potential described by Dalgic et al. 
to calculate the SSF of Cs and Rb in the realm of the 
liquid-plasma transition [37]. The author of [37J showed 
that calculated data were in agreement with the measured 
SSF. 

In Fig. S] the e — i and i — i HGK, screened HGK and 
i — i Dalgic et al. potentials are presented for comparison. 
One can easily see that with the growth of r ee (defined in 
the section |2~21 the curves shift in the direction of its low 
absolute values. We presume that this occurs due to the 
increasing role of screening effects. 



2.2 Static Structure Factors 

Within the framework of the density response formalism 
for two-component plasmas, we calculated the screened 
HGK interaction potentials described above using the semi- 
classical approach suggested by Arkhipov et al. [25]. Our 
approach is based on the HGK pseudopotential model for 



m r +m s 

reduced mass of the interacting pair r — s, and r , s = e 
(electrons) or i (ions). The effective temperature T' rs is 
given by , 

T , = m r T' s + m s T' r 



where T' e = (T|+T|)V 2 withT, = T F /(l. 3251-0.1779^), 
where r s = r a /r B , T F = fr 2 (37r 2 n e ) 2/3 / (2k B m e ) and TV = 
(T? + loTl) 1 / 2 , T D = Q pl h/k B , 7o = 0.152 is the Bohm- 
Staver relation for the Debye temperature with il 2 t = 

^pi/i^+^De/k 2 ), uj p i = y/ ze 2 n e / (eo m i) with n%i being the 
ion mass, k Be — W e 2 n e / '(e$k B T£) is the Debye wave num- 
ber for the electron fluid (T D w O.WeV, T F w 14.5eV for 
Be 2+ ). This definition of the effective temperature allows 
to extend the fluctuation-dissipation theorem to nonequi- 
librium systems and is the input value for the partial static 
structure factors S rs (k), the Fourier transform of the pair 
distribution functions h rs (r) — g r s(r) — 1 |29) : 



(12) S rs (k) = 5 rs + ^/n r n s h rs (k) ~ 5 r 



(13) 



Since, in the Debye model, the phonon modes with 
wavelengths up to a fraction of the lattice spacing are 
considered, in [T] it is set k = k max = (2j z) x l z kp with 
kp — (37r 2 n e ) 1 / 3 Fermi wave number. Due to the large 
mass difference between ions and electrons, T' ei = T' ee . All 
the parameters considered here are beyond the degenera- 
tion border (n e \ 3 e < 1). 

The partial static structure factors of the system are 
defined as the static (equal-time) correlation functions of 
the Fourier components of the microscopic partial charge 
densities 1351: 



S rs (k) = ^<p r (k)p s (-k)>, 



where N number of ions and 

N 



p r (k) = ^exp(ik- r. t 



(14) 



(15) 



i=l 



A linear combination of the partial structure factors which 
is of high importance, is the charge-charge structure fac- 
tor defined as 



S zz (k) 



1 



See (fc) - 2^S ei (k) + zSu(k) 



(16) 



where p z — p l (k) — p e (k). 

As described in [T] by Gregori et al., the fluctuation- 
dissipation theorem may still be a valid approximation 
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r ee = 0.606,4: Screened HGK at T e = T 4 = 13eV, T e ' = 18.65, T ee = 0.782, 5: HGK. The plasma parameters are described in 
the section [ 



even under nonequilibrium conditions if the temperature 
relaxation is slow compared to the electron density fluc- 
tuation time scale. A common condition in experimental 
plasmas for this to occur is when m t >> m e so that the 
coupling between the two-components takes place at suffi- 
ciently low frequencies. Using a two-component hypernet- 
ted chain (HNC) approximation scheme, Seuferling et al. 
[2"5] have shown that the static response under the condi- 
tions of the non-LTE takes the form: 

S rs (k) = 5 rs - ^^$ rs (k) - 5 er 6 es &- i)lS&. Sii {k) 

(17) 

where q(k) = y/zS e i(k) / Su(k) and for <P rs the expression 
(jTlfTUl) was used. 

In Figures[S](a), (b), (c), (d) the static structure factors 
S rs (k) for a beryllium plasma with the introduced above 
different temperatures T% = T e , T t = 0.5 • T e , T% = 0.2 • 
T e and the coupling parameters r ee — e 2 /(4eofcs7 1 e 'r ee ), 
r u = z 2 e 2 /{Ae k B Tlr u ) with r u = (3/4t™ 2 ) 1/3 , r ee = 
(3/47rn e ) 1 / 3 are shown. For typical conditions found in 
laser plasma experiments with solid density beryllium, we 
have n e m 2.5 • 10 23 cm~ 3 and z w 2. This gives T F sa 
lA.beV and Td ~ 0.l7eV. In Fig.[S](c) a minimum arises 
which defines the size of an ion core. Notice that in the 
Figure the minimum becomes less pronounced when the 
coupling increases. 

In a screened OCP the effective response of the medium 
is described by the charge-charge correlation function as 
given by Gregori et al. |12) : 



It is of high interest to study the influence of the ion 
structure on the static structure factors. For this reason 
in Figures [5] (a)-(d) and further we compare the SSF, ob- 
tained from equations ([T7]) and (|18D with the help of the 
screened HGK potentials ([7l fl"0|) . with the corresponding 
SSF obtained using the screened Deutsch potential |28j . 
on the basis of the TCP hypernetted-chain approximation 
developed for the case of non-LTE by Seuferling et al. [2T?] 
and further extended for DSF by Gregori et al. pQ. 

In Fig. [5] the static charge-charge structure factor ([15]) 
for a beryllium plasma with n e ~ 2.5 • 10 23 cm -3 , z ~ 2, 
T e = 20eV and T, = T e , T t = 0.5 • T e , Ti = 0.2 • T e 
is shown. In Figures [7] (a) - (d) we compare our results 
on the charge-charge SSF (TT5|) for alkali plasmas within 
the screened HGK model with the results obtained 
in the present work for alkali (hydrogen-like point charges 
(HLPC)) plasmas considered within the screened Deutsch 
model for various values of density and temperature. All 
curves obtained within the screened Deutsch model con- 
verge due to the negligible influence of an alkali ion mass 
on the wavelength A & scale entering the equations in [25] . 
As one can easily see with the growth of coupling the peaks 
become more pronounced and the difference among the 
curves becomes significant. We see that strong coupling 
and the onset of short-range order manifest themselves in 
S zz as a first localized peak, shown in an amplified scale, 
for different values of k! for every alkali species, and the 
position of the peaks shifts in the direction of small values 
of k' . This phenomenon was also reported in [12] . 
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Fig. 5. Static structure factors and screening charge q(k') for Be 2+ plasma at T e = 20eV, T' e = 24.0661/, z — 2 and n e = 
2.5 x 10 23 cm -3 . The set of filled symbols represents the screened Deutsch model obtained by Gregori et al.pQ , while the set 



of hollow symbols - the screened HGK model. Squares: Ti/T e 
r ee = 0.61). Triangles: T t /T e = 0.2 (P« = 11.57, T ee = 0.61). 



1 (r« = 2.31, r ee = 0.61). Circles: Ti/T e = 0.5 (r« = 4.63, 



3 The dynamic structure factor: the moment 
approach 

Extensive molecular-dynamic simulations of Coulomb sys- 
tems at a complete thermal equilibrium over a wide range 
of variation of the coupling parameter Fa [z — 1) and 
= Ep/ksT ( Ep is the Fermi energy) have been carried 
out by Hansen et al [35] • Hansen et al. studied dynamic 
and static properties of one-(OCP) and two-component 
plasmas and binary ionic mixtures. A new "moment ap- 
proach" based on exact relations and sum rules was sug- 
gested in [2] in order to calculate dynamic characteristics 
of OCP and of the charge-charge dynamic structure factor 
of model semiquantal two-component plasmas. This ap- 
proach proved to produce good agreement with the MD 
data of Hansen et al. Let 



JV 

E 



exp(tk-Ti r (t)). 



(19) 



be the Fourier components of the time dependent micro- 
scopic density of species r [35]. The corresponding dy- 
namic structure factors are the Fourier transforms of the 
density-density time correlation functions given by 



S rs (k,oj) 



1 



27TN 



<p r (k,t)p s (-k,0)>. (20) 



Alternatively, the charge-charge dynamic structure factor 
S zz {k,uj) can be defined via the fluctuation-dissipation 
theorem (FDT) [3] as 



S zz (k,uj) 



hlme 1 (fc,w) 



7r<£(fc)[l-exp(-/3fiw)]' 



(21) 



where <P(k) = e 2 /eok 2 and e _1 (fc, uS) is the inverse longitu- 
dinal dielectric function of the plasma. The charge-charge 
dynamic structure factor is directly related to the charge- 
charge static structure factor as follows : 
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Fig. 7. The charge-charge static structure factors S zz (|18p for alkali plasmas (Li + , Na + ,K + , Rb + , Cs + ) in a frame of the 
screened HGK model and results obtained in the present work for hydrogen-like plasmas in a frame of the screened Deutsch 
model on a base of Gregori et al.[T|. (a) T e = T; = 600007^, T' e = 60204A:, r ee = 0.398, Fa = 0.399; (b) T e = T t = 30000if, 
T' e = 30407^, r ee = 0.789, r u = 0.8; (c) T e = T t = 30000^, T' e = 31471K, F ee = 1.14, Fa = 1.2; (d) T e = T< = 30000_ft:, 
T' e = 37806if, F ee = 1.58, r la = 2. 



1 Z" 00 

"S^(fc) = ; / S zz (k,uj)duj 

n e + zni J_ 00 

_ S ee {k) - 2^S &i (k) + zS u (k) 



(22) 



where T' e =T' i =T e = T, T' ei = T' ee = T' e , n e = zn t (z = 1 
for hydrogen-like plasmas). 

In order to construct the inverse longitudinal dielectric 
function within the moment approach one needs to con- 
sider the frequency moments of the loss function 
— Ime (k,oj)/u>: 



C v {k) = -7T- 



L JTO£- X (fc,a;)dw, (23) 



Then the Nevanlinna formula of the classical theory of 
moments [40 expresses the response function 



e -1 (fc,w) = 1 + 



lo p 2 (uj + q) 



Uj(jjJ 2 — Ll>2 2 ) + l{^> 2 — ^l 2 ) ' 



(24) 



in terms of a Nevanlinna-class q — q(k, ui). The frequencies 
0J\ and uj2 are defined as respective ratios of the moments 
C u : 



ux 2 = C 2 /C = cj/[1 - e-^fc.O)]" 1 , 
u 2 2 = Ci/C 2 = u p 2 [l + Q(k)}, 



(25) 



where e~ l (k, 0) can be determined from the classical form 
(h — > 0) of the FDT (thermal equillibrium) eq. ([21~j) and 
the Kramers-Kronig relation [13) : 



Ree-\k,Lo) = 1 + -P.V. / K -^—^duj' 

7T 7_ oc W'-W 



(26) 
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Fig. 6. The charge-charge static structure factors S zz (|18ll 
for a beryllium plasma with n e w 2.5 • 10 23 cm -3 , z « 2, and 
T e = 2(W, T e ' = 24.06eV. The set of filled symbols represents 
the screened Deutsch model obtained here on a base of Gregori 
et al. [1] , while the set of hollow symbols - the screened HGK 
model. Squares: T;/T e = 1 (r tl = 2.31, r ee = 0.61). Circles: 
Ti/Te = 0.5 (r« = 4.63, T ee = 0.61). Triangles: T-/T e = 0.2 

(r« = 11.57, r ee = o.6i). 



In this way, we get the following expression : 

1. 2 



i?e£ _1 (fc,0) = 1 - 2S zz (k)- 



k 2 



(27) 



where flee -1 ^, 0) = e _1 ( fe >°) = £-*(&), 5^(fe) is denned 
by (|22[) . The function defining the second moment is given 
by 0: 

Q{k) = K(k) + L(k) + H. (28) 
It contains the kinetic contribution for a classical system: 



K(k) = 3(- 



(29) 



where kjj 2 = k^ 2 = n e e 2 / e^ksT . The contribution due 
to electron-ion HGK correlations is in our approach rep- 
resented by: 



H 



h el {r = 0) g ei (r = 0) - 1 



(30) 



Within the screened HGK model, the H in eq. (|30p 
can be approximated by —1/3 because we consider the ion 
structure, that means that the electron cannot approach 
the ion at r = distance. Note that g e i(r = 0) = is 
not exactly true, it is just a good approximation. The 
term L(k) takes into account the electronic correlations, 
in accordance to the approximation g e i(r = 0) = we 
calculated it for the Coulomb potential: 



L(k) 



1 



27T 2 7 



p 2 [S ee (p)-l]f(p,k)dp, (31) 



f(p,k) = 



12 



P 

4fc 2 



(fc 2 



8pk 3 



In 



p — k 



(32) 



In (|3Tj) the static structure factor is the one defined in 
(JTTJ) with the potentials given in (fT HTUj) . The authors of 
[3] suggested to approximate q(k,ui) by its static value 
q(k, 0) = ih(k), connected to the static value S zz (k, 0) of 
the dynamic structure factor through eq. (|2T1) : 



h(k) 



with 



7T^(k)uJ^S zz (k,0) 

S zz (k,0) ~ S°„(k,0), 



(33) 



where S°,(*,0) = ity^T^T 
dynamic factor takes the following form: 

s zz (k,Lj) = §h 

S zz (k,0) [l-exp(-^)] 

ajh 2 (k)uji 4 



(34) 

4T1 so that the normalized 



W 2 (W 2 -W 2 2 ) +/l 2 (fc)(w 2 -W! 2 ) 



with the more simplified expressions for h(k): 



h{k) 



e a y/2Trk B Tk 3 Lj p 2 (uj2 2 - ^i 2 ) 



, (35) 



(36) 



wi(fc), oj 2 (k): 



C2/C0 



2k D 2 S zz {k) ' 



u 2 2 = C4/C 2 - Wp 2 [l + tf(fc) + L(fc) - -], (37) 

In Figures [8] and [9] the DSF at a moderate tempera- 
ture T = 30000^ and concentrations n e = 1.741 • 10 20 - 
10 22 cm -3 but for the values of r ee used in [3] are shown. 
As one can see in Figures [5J |H1 the curves for alkali plas- 
mas are different from those given for the HLPC model 
[3], where the ion structure was not taken into account. 
In the case of alkali plasmas the curves split. This can 
be explained by that fact that alkali ion structure influ- 
ences the dynamic structure factor. In the Figures the 
position of the central peaks coincides but its intensity in 
alkali plasma is more pronounced. Positions of the plas- 
mon peaks are slightly shifted. We observe that the curves 
shift in the direction of lower values of k compared to the 
corresponding results of [3]- In Fig. [5] (b) the curves split 
into three very sharp peaks. Observe that with an increase 
of number of shell electrons from Li + to Cs + at low value 
of k the intensity of the lines grows while at higher value 
of k it diminishes. In Fig. [3] the position of the plasmon 
peaks shift in the direction of higher absolute value of ui, 
as compared to those in Fig. [8] This discrepancy could be 
also explained by that fact that the considered parame- 
ters are extreme, i. e., either high temperature or density 
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Fig. 8. Comparisons among the normalized dynamic structure factors of alkali plasmas (Li + , Na + , K + , Rb + , Cs + ) and the 
results of [3] for the HLPC model at k = 0.767/r ee , (a) T = 30000^, n e = 1.741 • 10 20 cm -3 , r ee = 0.5 and (b) T = 30000^, 
n e — 1.11 ■ 10 22 cm" 3 , r ee = 2. As the length scale we use the electron plasma frequency lj p — n e e 2 /eom e . 




and for r ee = 0.5, r s = 0.4 the degeneration condition is 
rie^ee = 0.335, while for r ee = 2, r s — 1 the degeneration 
condition is n e X ee = 0.678. 



4 Conclusions 

The e — e , e — i, i — i and charge-charge static structure 
factors have been calculated for alkali and Be 2+ plasmas 
using the method described and discussed by Gregori et 
al. in pp. The dynamic structure factors for alkali plas- 
mas have been calculated using the moment approach 2 , 
|3J. The screened Hellmann-Gurskii-Krasko potential, ob- 
tained on the basis of Bogolyubov's method, has been used 
taking into account not only the quantum-mechanical ef- 
fects but also the ion structure pi]. The results obtained 



within the screened HGK model for the static and dy- 
namic structure factors have been compared with those 
obtained by Gregori et al. and in [3] for the hydrogen-like 
point charges model. We have found small deviations (in 
the values of the SSFs) from results obtained by Gregori 
et al while there are significant differences between DSFs. 
Nevertheless, we have noticed that the present results are 
in a reasonable agreement with those of [3] at higher val- 
ues of k and with increasing k the curves damp while at 
lower values of k we observe sharp peaks also reported in 
[3]. At lower r ee the curves for Li + , Na + , K + , Rb + and 
Cs + do not differ while at higher r ee the curves split. As 
the number of shell electrons increases from Li + to Cs + 
at low k the intensity of the lines grows, while at higher 
k it diminishes. The difference is due to the short range 
forces which we took into account by the HGK model in 



12 



S. P. Sadykova et al.: Stat, and Dyn. Struc. Factors with Account of the Ion Structure 



comparison with the hydrogen-like point charges model. 
One should also take into account that we employed dif- 
ferent plasma parameters because at the parameters used 
in [3] the alkali plasmas with closed shell cannot exist. 
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